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Abstract 

I discuss a recent application of homotopy perturbation and Adomian decomposition 
methods to the linear and nonlinear Schrodinger equations. I propose a generaliza- 
tion of the procedure for the treatment of a wider class of problems. 



1 Introduction 

The Homotopy Perturbation Method (HPM) [1-15] (the reader should not 
mistake this abbreviation with that for the Hypervirial Perturbation Method 
[16]) is supposed to solve any linear and nonlinear problem in any field of 
theoretical physics. In a couple of papers I have shown that some of those 
results are useless and worthless [17,18]. My criticisms have not been welcomed 
by the referees (most probably homotopy devotees) and therefore they remain 
unpublished. 
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The purpose of this communication is to analyze a recent application of 
HPM and Adomian decomposition method (ADM) to the linear and nonlinear 
Schrodinger equations [15]. Since those equations are known to be relevant to 
many fields of physics, any approach to solve them is always welcome. 

2 Linear Schrodinger equation 

Sadighi and Ganji [15] consider the linear Schrodinger equation 

u t (x,t) + iu xx (x,t) = (1) 



where i 2 = — 1 and the subscripts indicate derivatives with respect to the vari- 
ables. We immediately recognize this differential equation as a dimensionless 
version of the Schrodinger equation for the free particle. Its propagator is well 
known and there is no reason for any special method to solve it. However, 
Sadighi and Ganji [15] think otherwise and apply the HPM to it. 

Most textbooks of quantum mechanics solve equation ([T]) for a physical square- 
integrable initial function u(x, 0) 

oo 

J \u{x,0)\ 2 dx < oo (2) 

— oo 

however an homotopy devotee does not care about such a small detail. In 
fact, Sadighi and Ganji [15] cleverly choose some particular unphysical initial 
functions, so that they solve the resulting equations much more easily. 

Any undergraduate student in a first course of quantum mechanics knows that 
if one chooses an unphysical initial function of the form 

u(x,0) = J2 c 3 e ° jX ( 3 ) 
5 



2 



then the resulting unphysical solution will be 



u[ x 



a-jX—iott 

e 3 •? 



(4) 



Notice that the initial function (j3J) does not statisfy (J2J) 



The first example chosen by Sadighi and Ganji [15] 
u(x,0) = 1 + 2cosh(2x) 



(5) 



is a particular case of ([3]). If you think that it is a nonsensical quantum- 
mechanical state then you are not an homotopy devotee. It follows from equa- 
tion (H that 



u(x, t) — 1 + 2 cosh(2a;)e" 



-Ait 



(6) 



After a masterful application of the HPM Sadighi and Ganji [15] obtain 



u(x, t) = 1 + 2 cosh(2a;) 



1-41* + !^ + ... 



(7) 



and brilliantly conclude that the limit of this series is equation (jHJ). 

Example 2 is even simpler than example 1; after some algebraic manipulation 
Sadighi and Ganji [15] are led to the revelation that when 



u(x, 0) = e' 



?A 1 



(8) 



then 



u(x, t) 



l + 9it+ + ... 



3 3(x+3ii) 



(9) 



which is obviously a particular case of equation (j3J). I insist that if you think 
that this solution is completely useless from a physical point of view, then you 
are not an homotopy devotee. 
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3 Nonlinear Schrodinger equation 



The nonlinear Schrodinger equation 

iu t + u xx + -y\u\ 2 u = (10) 



is a much more interesting and challenging problem. However, Sadighi and 
Ganji [15] cleverly restrict themselves to unphysical solutions with unit mod- 
ulus \u\ = 1 so that the problem reduces to a trivial linear differential equation 
with constant coefficients: 

iut + u xx + yu = (11) 



If you think that this restriction is ridiculous then I tell you that you are not 
an homotopy devotee. 

Any student in a first course of calculus will immediately realize that a par- 
ticular solution is 

u(x,t) = e iax e i{ ^- a2)t (12) 

When a = 1 and 7 = 2 we obtain 

u(x,t) = e l{x+t) (13) 



which is exactly the example 3 of Sadighi and Ganji [15]. After masterfully 
solving the tedious HPM equations Sadighi and Ganji [15] obtain 

(it) 2 



u(x, t) 



it 



2! 



+ . . . 



(14) 



which obviously converges to (TIBl . 
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When a = 1 and 7 = 

u(x, t) = e i(x - 3t) 



—2 equation (fT2|) reduces to the solution of example 4: 

(15) 



As the reader may have guessed, after solving the HPM equations Sadighi and 
Ganji [15] obtain the power series 

(3it) 2 



u(x,t) = e lx 



1-3U + 



2! 



(16) 



that clearly converges towards (ITo|) . 

It is pointless telling an homotopy devotee that all those calculations are 
useless and worthless, he is content with showing that his calculation yields 
the exact result known to everybody. 

The reader should not think that Sadighi and Ganji [15] merely restrict to the 
calculation of the power series of exponential functions e lbt by means of HPM. 
They also do exactly the same by means of another of the most fashionable 
approaches: the Adomian decomposition method (ADM) [15]. 



4 Conclusions 



The discussion of the preceding sections show that Sadighi and Ganji [15] 
study the simplest linear and nonlinear Schrodinger equations. They choose 
unphysical solutions that do not correspond to actual quantum-mechanical 
states but are more tractable. One such choice converts the nonlinear Schrodinger 
equation into a linear one and consequently the authors never solve the nonlin- 
ear problem. Then they apply two approximate methods, the HPM and ADM, 
and obtain a power series approximation to the solutions. As shown above, 
the calculation reduces to the Taylor expansion of exponential functions of the 
form e M about t = 0. If the reader thinks that it is not a great achievement I 
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will prove him/her wrong. First of all, notice that the authors do not use the 
obvious formula 

(17) 

t=o 

but apply two cumbersome methods in order to obtain this expansion. Second, 
they managed to publish their remarkably useless and nonsensical contribution 
in a research journal [15]. In another contribution I will analyze how one can 
do it. 

Finally, I will like to propose a generalization of the great achievement of 
Sadighi and Ganji [15]. Instead of the simple linear equation ([T]) one may 
consider the more general one 

u t + iAu = (18) 



u(x, t) 



v V d j u 



where A is a linear operator. According to the textbooks on quantum mechan- 
ics, its formal solution is 

u = e - itA u (19) 



where uq = u(t = 0). A successful application of the HPM or ADM will lead 
to the series 



u 



i (-itA 2 

1-UA+ - r-t- + ... 

2! 



(20) 



If you do not succeed in publishing your results you may try a simpler problem; 
for example, when 



Au = CLUq 



(21) 
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you have 



u = e 



-ita 



U 



(22) 



and 



\ . (-ita) 2 
u= l-ita+ v - ' + 



U 



(23) 



Hopefully, the referee (most probably an homotopy devotee) will understand 
your simplified problem and then accept your manuscript. 

If your manuscript is accepted then you become a member of the homotopy 
club and then you can try something bolder, for example 

uo = °kfk (24) 
k 

where 

Afk = a k f k (25) 



I leave the result as an exercise for the future homotopy devotee. 

Notice that you can thus write a paper for every solvable or unsolvable problem 
which you can think of. Most probably, all of them will be accepted in any one 
of the homotopy journals [1-15] (unless you dare to criticize another homotopy 
result). 
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